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Abstract 

We consider a special case of the n-component cubic model on the square lattice, for which an 
expansion exists in Ising-like graphs. We construct a transfer matrix and perform a finite-size- 
scaling analysis to determine the critical points for several values of n. Furthermore we determine 
several universal quantities, including three critical exponents. For n < 2, these results agree well 
with the theoretical predictions for the critical 0(n) branch. This model is also a special case of 
the {Na,Np) model of Domany and Riedel. It appears that the self-dual plane of the latter model 
contains the exactly known critical points of the n = 1 and 2 cubic models. For this reason we 
have checked whether this is also the case for 1 < n < 2. However, this possibility is excluded by 
our numerical results. 
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I. INTRODUCTION 



The n-component cubic model can be defined in terms of vector spins tliat are restricted 
to lie along one of n Cartesian axes, but are free to assume the positive or negative direction. 
Only one vector component is nonzero; it is normalized to be ±1. The model can be described 
by the reduced Hamiltonian 

n/k^T = - Y.[Ks, ■ Sj + M{si ■ Sj)^] , (1) 

where the index i of the cubic spin Si refers to the sites of the square lattice. The sum on 
{ij) stands for all nearest-neighbor pairs. This model obviously combines Potts degrees of 
freedom (the choice of the Cartesian axis, which is subject to permutation symmetry) with 
Ising degrees of freedom which specify the sign of the nonzero component. This is more 
explicit in the following form of the Hamiltonian 

H/k^T = - Y.iK^^^3 + M)5^.r, , (2) 

where we represent the sign of the nonzero component of §1 hj Si = ±1 and its Cartesian 
axis number by = 1, 2, ■ ■ ■ , n. The corresponding bond weight Wij can be rewritten as 

Wij = exp[{KSiSj + M)6r^rj] 

= 1 + 6r,r,{exp{KsiSj + M) - 1} 

= 1 + 6r^r,n{v + XSiSj} 

= j^iSr^rM^ + XS.Sjjp. (3) 

In the third line we have used the definitions 

e*^ cosh K ~ 1 sinh K , . 

V = and X = , (4) 

n n 

and in the fourth line we have introduced bond variables bij = or 1, and the summand is 

subject to the rule 0° = 1. Thus the partition sum assumes the form 

1 

Zcuh = E n E l^nrM^ + ^■'^iSj}]'''' ■ (5) 

Application of the Kasteleyn-Fortuin mapping jlj involves execution of the sum on the 
Potts-like variables {r}. This leads to 

^cub = E [ n E + xs,sj}''^^]n^^ , (6) 

{s} {ij) fe.,=0 



where is the number of sites of the lattice. Note that each bond (by which we mean a 
bond variable bij = 1) contributes, through the Kronecker 6, also a factor 1/n, unless it 
closes a loop. The latter condition is accounted for by the factor n"', where ni is the number 
of loops formed by the bond variables. Each configuration of bond variables bij defines a 
graph on the square lattice covering those and only those edges for which b^j = 1. For the 
special case 

coshK = e-^^ or v = 0, (7) 
the partition sum reduces, after execution of the sum on {s}, to 

Zcub = (2n)^^a;"''n"', (8) 

{b} 

where the sum on {b} contains only 'even' graphs in which every site is connected to an 
even number of bonds bij = 1. The odd graphs, while included in the sums on the bij in 
Eq. do not survive the sum on {s}. The number of bonds in the graph {b} is denoted 
nb = J2{ij) bij- 

The cubic loop model described by Eq. (jHl), is subject to the restriction \xn\ < 1 because 
of Eq. (0). This model, and the model defined by Eq. with M = 0, (both on the square 
lattice) have already been investigated by a transfer- matrix method . That work included 
a determination of the temperature exponent. The results for n < 2 were in agreement 

n n 

with a conjecture of Cardy and Hamber [3|| and the analysis of Nienhuis jj]. However, the 
latter work applies to the n-component cubic model on the honeycomb lattice, of which 
the graph expansion reduces to that of the 0(n) model. For other lattices, no such direct 
correspondence exists, and the relevance of the results of Ref. I4 for the cubic model on the 
square lattice thus needs further justification. This is provided by renormalization arguments 
that predict that cubic anisotropy is irrelevant [5j for n < 2. However, for n = 2 it is marginal 



and indeed the temperature exponent does not agree with the 0(2) model 

The work of Ref. \^ did not include results for the interval 1 < n < 2, where the emphasis 
lies of the present work. Furthermore, the present analysis will purportedly yield more 
information about universal quantities. 

The outline of this paper is as follows. Section |H] discusses a duality transformation in a 
three-dimensional parameter space containing the cubic model Eq. (jH)). This transformation 
suggests a possible form of the critical line of the cubic model. However, transfer-matrix 
calculations, defined in Sec. IIIIl yield numerical results listed in Sec. IIVI which show that 



this form is not applicable. The universal properties of the cubic model are investigated in 
Sec. El The paper concludes with a short discussion in Sec. IVII 



II. SELF-DUAL PLANE OF THE (iVa,A^/3) MODEL 



The partition function of the (Na^Nfs) model [6[ is defined as 

m n 
a=l T=l edges 

where m = N^, n = N^, and 

B{(T, T\a', t') = eXp[KoSaa'SrT' + Ki{l - Saa')SrT' + 
K26„a'{l - Srr') + ^3(1 - ' Srr')] • 



(9) 



(10) 



Just as the cubic model, the model can be viewed as having two Potts spins a and r on each 
lattice site, with allowed values a = 1, 2, ■ ■ ■ , m and r = 1, 2, ■ ■ ■ , n. These spins interact 
with nearest-neighbor couplings according to the Boltzmann weights (fTIUl . Note that these 
weights are elements of the (mn) x (mn) matrix 



B = e^' 



(11) 



where /„ and J„ are n x n matrices 
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(12) 



We consider a model on a planar lattice with open boundary conditions. The partition 
sum possesses the duality relation Q 



(--)^-^^<l)(e^«,e^Se^J,e^3 



(13) 



where N£) is the number of sites of the dual lattice, and Z^^'^ the dual partition function, 
with Boltzmann weights also given by Eq. ([111)) but with the Ki replaced by the dual weights 



K*. The latter weights were shown W to be equal to the eigenvalues of the matrix B. For 
the (m, n) model, these are 



el 



ye ^ J 



( 



1 m — 1 n—1 (m — — 1) 
1 -1 n-1 
1 m — 1 —1 — (m — 1) 
1-1-1 1 



(14) 



ye 



Note that p4|) implies e^^**-** = {mn)e^\ and hence (Z*^^^)*^-^) = Z, because the number of 
lattice edges E satisfies Euler's relation E = N + Nd — 2 and each lattice edge contributes 
a factor mn. Eq. (jl3p then shows that the powers of mn cancel after a pair of duality 
transformations. 
With the notation 



1,2,3) 



(15) 



the Boltzmann factor can be written as 



B(a,T\a',r') 



e "X 



X' 



(16) 



Phase transitions will naturally occur in the x = {xi,X2,X3) parameter space. Thus, the 
free energy of the (m, n) model is expected to be nonanalytic on a corresponding surface S. 
Generally, as the temperature variable in a given model varies, the point {xi,X2,Xs) traces 
out a certain "thermodynamic" path F in x, and the model exhibits a transition whenever 
F crosses S. 

Using (fT^ we find the transformation 

1 



X, 



X n 



[1 - Xi + {n - l)x2 - {n- 1)0:3] 
[1 — (m — l)a;i — X2 — (m — 1)0:3] 



A 

1 

A 

— [1 ~Xi-X2 + X3] 



where 



A = 1 + (m — l)xi + (n — 1)0:2 + (m — 1)(?7. — 1)0:3 



(17) 



(18) 



The square lattice maps onto itself under the dual transformation, so that the free energies 
at two points Xi and x* satisfying Eq. (fT7|) are related. One also verifies that the subspace 



A = ^/rrm (19) 

is self-dual under the transformation Eq. ()17|) (but a point in this plane does, in general, 
not map on itself). 

We note that the partition sum of the cubic loop model described by Eq. (0) with f = 0, 
i.e., the model of Eq. (jSJ, can be written as 

Zcnh= W[l+nXSiSj5r,rj\- (20) 

{s},{r} (ij) 

By introducing ak = {sk + 3)/2, and the identities 

S^SJ = 26^^^^ - 1 (21) 

we find the equivalence 

Zcnh = ^(2,n)(e^°, e^S e^^) (22) 

with 

= l + nx, e^^ = l-nx, e^^ = e^^ = 1 (23) 

or 

1 — nx 1 1 . ^. 

= -r, ' ^2 = — , X3 = — . 24 

1 + nx 1 + nx 1 + nx 

This specifies the mapping of the n-component cubic loop model on the (A^^ = 2, N/s = n) 
model. Using the dual transformation (fT7|) . this gives rise to 

x^ = x^ = X, X2 ~~ . (2^) 

Hence the dual thermodynamic path F of the cubic loop model is the straight line connecting 
(0, 0, 0) and (1, 0, 1), a result valid for all n. 

For two special cases, namely n = 1 and n = 2, the critical point of the cubic loop model 
sits at the intersection of the critical surface S and the thermodynamic path F in the X2 = 
plane. 

We first consider the case n = 1, or (5^-^^-^ = 1, in which the model simply reduces to the 
Ising model. The Boltzmann factor (fTB|) assumes the form 

B{a,T\a\T')=e''^x\~^-^' (26) 



for any values of X2,X3. The critical surface is thus Xi = \/2 — 1. This plane intersects the 
thermodynamic path T at xi = X3 = x = 1, so that the critical point of the cubic loop 
model is Xc = \pl — 1. 

For the case n = 2, the model, i.e., the (2,2) model, is the well-known Ashkin- Teller 
model ^. The shape and location of S have been discussed by Wu and Lin P]. The 
thermodynamic path F crosses the critical surface S at xi = X3 = 1/2 within the X2 = 
plane p!]. Thus the critical point of the cubic loop model occurs at Xc = 1/2. 

While it is known that most of the self-dual plane of the {Na, Nf^) model is noncritical, it 
is interesting that for both n = 1 and n = 2, the critical points of the cubic loop model are 
actually located in the self-dual plane. If this is true for general 1 < n < 2, we would have 



l + xi + {n- l)x3 = V2n (27) 

on the X2 = plane, and Xi = X3 = x, because the critical point lies on the thermodynamic 
path F. Thus the critical value of x would be 



x^{n) = {V2n -l)/n for 1 < n < 2 . (28) 
This possibility will be investigated numerically in Sec. IIVI 

III. THE TRANSFER MATRIX 



The transfer-matrix method used here is related to that used in Ref. |2, and it uses in 
addition some of the techniques described in Refs. 10 and 3 for the random-cluster and the 
0{n) model respectively. The full description of the transfer matrix is somewhat elaborate, 
and here we only provide a general outline, supplemented with more detailed information 
where the procedure is different from those in the references given. 



As in Ref. |2, the transfer matrix is constructed on the basis of a graph representation 
of the cubic model that allows n to be non-integer. However, the present work is restricted 
to the case f = 0, so that the graphs are restricted to be even. This allows the use, given 
a system size, of a smaller transfer matrix than that used in Ref. 2I We define the model 
on an L X m lattice Cm wrapped on a cylinder, such that the finite-size parameter L is the 
circumference of the cylinder. The definition of the transfer matrix can be illustrated by 
appending row m + 1 and determining how this affects the partition sum of the model. The 



lattice Cm has an open end at row m; there are L 'danghng' edges that will serve to connect 
to row m + 1 later. Whereas the partition sum Eq. © allows only for closed loops, the 
bond configuration on may be considered a part of a larger graph so that we allow the 
dangling edges to be occupied by loop segments. But all sites are still restricted to connect 
to an even number of bonds. For the construction of the transfer matrix, we need a coding 
of all possible ways (called connectivities) that the dangling bonds can be connected by the 
graph {6} on Cm- This coding assigns a unique integer 1,2,---, which will serve as the 
transfer-matrix index, to each connectivity. Some of the dangling edges may be empty, i.e., 
hij = 0. The remaining dangling edges, i.e., the dangling bonds, form a 'dense connectivity' 
without vacant positions. Note that these dense connectivities satisfy a 'well-nestedness' 
principle which asserts that, if positions i and k are connected, and j and / are connected, 
with i < j < k < I, then all i,j,k,l must be connected. Thus, these dense connectivities 
form a subset of the 'random cluster' or 'Whitney' connectivities defined in Ref. The 
number of dangling bonds that are connected by a path of bonds must always be even for 
the cubic model; this restriction does not apply to those of Ref. Q. By simply excluding 
the odd connectivities, we thus find coding and decoding algorithms for the dense cubic 
connectivities. The coding of the general cubic connectivities including vacant positions 
then follows analogous to the case of the 'magnetic' connectivities in Ref. [l^ 

This coding allows one to divide the partition sum of the m-row system in contri- 
butions corresponding with different connectivities: 

^(-) _ ^ ^(m) _ (29) 

a 

Next, we append a new row Im+i to the lattice: Cm+i = 'Cm^lm+i, and express the restricted 
partition sums as a linear combination of the zj^^ . This is possible because the weight 

due to the newly appended row is completely determined by the bond variables connecting 
to the appended vertices and the 'old' connectivity (3, and this information also determines 
the 'new' connectivity a. We use &m+i to denote the 2L appended bond variables, and 
&m+i) to denote the function that determines a. The weight factor associated with the 
new row satisfies 

t/;(/3,6„+i) = (2n)^a;^"^n^"S (30) 
where Arib is the number of appended bonds and Ani is the number of loops closed by these 



bonds. The recursion connecting the restricted sums is 

^(M+l)^^y^^^W^ (31) 

in which the transfer matrix Tap is defined by 

Tap = {2nf • (32) 

bm + l 

In actual calculations, the transfer matrix is represented as the product of L sparse matrices, 
each of which appends one new vertex of the (m + l)-th row. The first vertex of a new row 
increases the number of dangling edges to L + 2, so that the sparse matrices assume a larger 
size than Tap- After appending the last vertex of that row, the number of dangling edges 
decreases to L. This technical point was described in some detail for the related case of the 
0{n) model on the square lattice 

The sparse-matrix technique makes it unnecessary to store the full transfer matrix T. 
Some of its eigenvalues can be obtained by repeated multiplication of a vector by T, and 
a„a,y=. of the .exulting vector sequence. Stnce . not general, we nsed 

the method of projection to a Hessenberg matrix as described in Ref. UJJ. We restricted the 
calculations to vectors with 'translation symmetry', i.e., vectors that are invariant under a 
permutation of connectivities corresponding with a cyclic permutation of the dangling edges. 
In general, the largest eigenvalue A^^ determines the free energy f{L) per site in the limit 
of an infinitely long cylinder (m oo): 

/(L) = L-MnAf . (33) 

Furthermore, the next largest eigenvalues A^'' (z = 1, 2, ■ ■ ■) determine the correlation lengths 
^j(L) of various types of correlation functions. The latter types depend on the symmetry 
of the corresponding eigenvector and on possible modifications of T. In particular, the 
correlation length ^t{L) of the energy-energy correlation function is determined by the gap 
between the two largest eigenvalues: 

C^(L)=ln(Ai°VA?V (34) 

For the cubic model, magnetic correlations can be represented, in analogy with the 0(?7,) 
model, by graphs with odd vertices on the correlated sites. For the present model that 
means sites connected to one or three bonds. The two correlated sites, which are placed far 



apart in the length direction of the cyhnder, must be connected by the graph, i.e., belong 
to the same component of the graph. This additional component does not follow the rules 
of 'evenness' listed earlier. The number of danghng bonds at the open end of the cylinder 
(between the correlated sites) connecting to the additional odd vertex must be odd. To 
describe such 'magnetic' graphs we define a new set of connectivities in which one group 
of connected dangling bonds is odd. This leads to a modified transfer matrix, which may 
alternatively be interpreted as the magnetic sector of a larger transfer matrix whose vector 
space includes both even and odd connectivities. The gap between A^'* and the largest 
eigenvalue A^-* in the magnetic sector determines the magnetic correlation length ^^(L) as 



A different type of magnetic gap is associated with the density of the loops spanning the 
circumference of the cylinder. The weight of these loops is modified by assigning a bond 
weight —X to one bond in each row. A loop spanning the circumference contains an odd 
number of these modified bonds, and its weight thus changes sign. All other loops contain 
an even number of modified bonds and their weights are thus unchanged. We denote the 
largest eigenvalue of the modified transfer matrix by A^p . It determines the length scale 
^m{L) that may be associated with the effect of an 'antiferromagnetic seam' running along 
the cylinder. For the critical Ising case n — 1, both magnetic length scales are related by 
duality, but this is not so for general n. The corresponding length scale is given by 



In the actual transfer- matrix calculations, we have used finite sizes up to L = 15 in the 
nonmagnetic sector, which then has dimensionality 2 004032, and up to L = 14 in the 
magnetic sector, which then has dimensionality 3 856582. 

IV. DETERMINATION OF THE CRITICAL LINE OF THE CUBIC MODEL 

The asymptotic behavior of the magnetic correlation length ^h{L) near a critical point 
can be expressed in terms of the scaled gap 



C^(L)=ln(AS°VA?^)- 



(35) 




(36) 



Xh{t,u,L) 



L 



(37) 



where t parametrizes the distance to the critical point, and u represents an irrelevant field. 
Renormalization arguments 12 1, scaling and conformal invariance predict that for 
large L 

Xh{t, u, L)c^Xh + aiUH + biLy-u + ■ ■ ■ , (38) 

where is the magnetic scaling dimension, yt the temperature exponent, and Uu the ex- 
ponent of the field m, and ai and hi are unknown amplitudes. Further corrections may also 
be present. Since we have an algorithm available that calculates Xh{t,u,L) (with t and u 
expressed as a function of x), we can estimate the critical point by numerically solving x in 
the equation 

X,,{x,L)=Xf,{x,L + l), (39) 

which is a form of 'phenomenological renormalization' After substitution of Eq. (j38|) 
one finds that, at the solution, t and u satisfy 

t oc uL^^-y' . (40) 

Since yt > and ?/„ < 0, we expect that t — > for L ^ oo, i.e., the solutions of Eq. (j39|) . 
which we denote x^^\L), converge to the critical point. These solutions were fitted by solving 
for x^^\L), c^^\L) and yu — yt in the three following equations with L' = L, L — 1, and L — 2: 

(L') = x(^) (L) + c(^) {L)L'y--y' , (41) 

which leads to a sequence x'^^^L) that is shorter than the original sequence x^^\L) but 
usually shows faster apparent convergence. Another iteration step can be attempted on the 
basis of 

x(^) {L') = (L) + c(') {L)L'y'--y' , (42) 

which may lead to even better estimates of the critical point. 

A similar analysis of the critical point can be performed on the basis of the scaled interface 
gap 

Xmit,u,L) = -—^ n' (^^) 

using the same type of fits as for the scaled magnetic gap. 

We have also attempted to find solutions of Eq. with C,h replaced by the energy- 
energy correlation length, but here complications arise. The functions Xt{x,L) typically 
display an extremum near the critical point, and solutions of the scaling equation Eq. ()39p . 
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FIG. 1: Critical points of the v = cubic model. The data points show the numerical results of 
the transfer-matrix analysis, and the curve is added as a guide to the eye. The estimated error 
bars are smaller than the size of the symbols. 

with Xt instead of X^, do not always exist. In particular for n > 1 we did not obtain a 
satisfactory set of solutions, and we have not pursued this way to obtain further data for the 
critical point. Instead, we located the extremum of Xt{x,L) as a function of x. The finite- 
size-scaling equation for the correlation length indicates that this extremum will converge 
to the critical point. 

The estimated critical points are shown in Fig. ^ For 1 < n < 2 they do not agree with 
Eq. ()28|) . For instance, that equation would predict Xc = 0.48803 • ■ ■ for n = 1.5, which is 
incompatible with the numerical result (see also Table H}. Thus we have to conclude that 
Eq. ()28|) . which does indeed lack a solid basis, is not valid for all values of n. 

Several modified fitting procedures were applied. Assuming that the cubic model reduces 
to the 0{n) universality class, we have analytic evidence for the values of yu and y'^ as a 
function of n. First, according to the renormalization scenario, the cubicperturbation of 
the 0(n) symmetry corresponds with an irrelevant field with an exponent jsf 



where cosing) = —n/2 and I < g < 2. Corrections due to this field are expected to be 
serious for n —>■ 2, since yc 0. Second, irrelevant temperature-like fields may correspond 
with some scalar operators whose dimensions are entries Xi^q with g > 3 in the Kac table 




(44) 



Vc = 



2g 



1 o 



TABLE I: Critical points x^, as determined from the scaling formulas for the magnetic and interface 
correlation length for system sizes L and L + 1, and from the extrema of the energy-like correlation 
length as x is varied. The estimated numerical uncertainty in the last decimal place is shown in 
parentheses. The 'best estimates' are based on the results in the two preceding columns and on 
an analysis of the minima in the functions Xt{x,L) as described in the text. For n = 1 we find 
accurate agreement with the exact result Xc = \/2 — 1, and for n = 2 with Xc = 1/2. 
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with m = l/(g — 1) for the cubic model [l8|. For g = 5 we find an irrelevant exponent 

6(7-12 



9 



(45) 



(46) 



which has small values for n > but becomes marginal when n —2. The final results, 
and their estimated errors, were obtained from the analyses of the three types of scaled gap, 
and the degree of consistency between different types of fits mentioned above. The best 
estimates obtained from Xh{t,u,L) and Xm{t,u, L), and the overall best estimates which 
also include data from Xt{t, u, L), are shown in Table |l] 



V. UNIVERSAL PROPERTIES OF THE CUBIC MODEL 



The asymptotic finite-size dependence of tlie free energy per site at the critical point is 



m 



19| 



/(L) ^ /(oo) + ^ , (47) 

where c is the conformal anomaly of the model, which characterizes universality classes and 
determines sets of critical exponents 



2l|. We have calculated the finite size data for 
the free energy at the extrapolated critical points, and estimated c as c^^\L) from the free 
energy density for two consecutive system sizes by solving 

c«(L) = 6[/(L) - f{L + WW/L' - 1/{L + 1)^] . (48) 

This leads to a sequence of estimates of c that can be extrapolated by means of power-law fits, 
analogous to the procedure used to determine the critical points. After a second iteration 
step, the estimates of c seem almost converged, in the sense that the results for the largest 
few values of L display differences of only a few times 10^^. But at the same time these 
data display a shallow extremum (except for n = 1 and 2, where the apparent convergence 
is much better), so that it is difficult to estimate the uncertainty in the extrapolated results. 
Because the first iteration step shows that the finite-size dependence of the estimates of c is 
close to —2, we have also applied iteration steps with the exponent fixed at this value. The 
results were similar to those obtained with free exponents, and again displayed a shallow 
extremum. Under these circumstances we made a crude error estimate of ten times the 
difference between the two estimates for the largest available L values, after two iteration 
steps. The best estimate was taken by extrapolating the last two estimates, using again ten 
times the aforementioned difference. A better apparent convergence was found when a fixed 
exponent ?/c ~ 2 was used in the second iteration step. The results are shown in Table ITTl 
The numerical errors were estimated from the finite-size dependence of the results of the 
last fit procedure, except for n = 1.9, where the error bars of both fit procedures did not 
overlap and we used the difference between both types of fit instead. Most of our results are 
in good agreement with the theoretical values: 



9 



(49) 

which follow after the substitution of the formula m = l/{g — 1) in the relation [2^ 
between m and c, i.e., c = 1 — 6/[m(m + 1)]. But the result for n = 1.9 does not agree well 



TABLE II: Conformal anomaly and temperature scaling dimension Xf as determined by transfer- 
matrix calculations described in the text. Estimated error margins in the last decimal place are 
given in parentheses. For comparison, we include the theoretical values of the 0(n) model for 
n < 2, and of the Ashkin- Teller model for n = 2. The numerical results are indicated by 'num', 
theoretical values by 'theory'. 
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0.7418 (1) 
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1.251 (2) 


1.25189 
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0.7912 (1) 
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1.319 (2) 


1.31996 
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0.8410 (1) 


0.84096 


1.405 (5) 


1.39962 
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0.8920 (2) 
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1.49 (2) 


1.49783 
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0.9464 (9) 


0.94432 


1.60 (5) 


1.63279 


2.0 


1.00000 (1) 


1 


1.500000 (1) 
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with the theoretical value; we note that the small value of the cubic crossover exponent, 
which becomes marginal at n = 2, may well lead to imprecise results and error estimates. 

Next, we analyze the results for the magnetic gaps. After substitution of the solutions of 
Eq. (jSni), which behave as Eq. (jinjl . in Eq. (p!Hjl . one finds that the magnetic scaled gaps at 
the solutions converge to Xh as: 

XhiL) =Xh + ruiy- + ■ ■ ■ , (50) 

where r is an unknown constant. The magnetic scaled gaps at the solutions of the scaling 
equation in the previous section were already available. They were fitted using a similar 
procedure as used for the determination of the critical points. We found that the leading 
irrelevant exponent was consistent with the predicted cubic perturbation exponent given in 
Eq. ()44|) . and we accordingly treated tju as a known parameter in the fits. The extrapolated 
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TABLE III: The magnetic dimension and the interface dimension as extrapolated from 
their values at the solutions of the scaling equation of the correlation length. Estimated numerical 
uncertainty in the last decimal place are given in parentheses. For comparison, we include the 
theoretical values of the 0(n) model for n < 2, and of the Ashkin- Teller model for n = 2. The 
numerical results are indicated by '(num)', the theoretical values by '(theory)'. 
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n 194QQ (^) 
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1.1 


0.12647 (1) 


0.12668 


0.14105 (5) 


0.14101 


1.2 


0.1280 (1) 


0.12826 


0.15814 (5) 


0.15815 


1.3 


0.1296 (1) 


0.12973 


0.1767 (1) 


0.17668 


1.4 


0.1307 (1) 


0.13107 


0.1969 (2) 


0.19695 


1.5 


0.1316 (1) 


0.13224 


0.2195 (5) 


0.21946 


1.6 


0.1320 (1) 


0.13319 


0.244 (1) 


0.24499 


1.7 


0.1335 (3) 


0.13382 


0.274 (2) 


0.27486 


1.8 


0.1308 (2) 


0.13393 


0.313 (2) 


0.31169 


1.9 


0.128 (1) 


0.13300 


0.350 (5) 


0.36230 


2.0 


0.12498 (1) 


1/8 


0.37500 (1) 


3/8 



magnetic scaling dimensions are shown in Table IIIIl Again, different fit procedures, with 
the finite-size exponent left free, and more iteration steps, were applied. The error estimates 
are based on the apparent convergence and on the differences between the various types of 
fit. The final results appear to agree with the theoretical values for the 0(n) universality 
class Q: 

X. = ^-il^. (51) 
where we note that the Ashkin- Teller model (or (2, 2) model) has the same magnetic scaling 



dimension 



22 



23| X;, = 1/8 as the 0(2) model. 
A similar analysis was performed on the scaled interface gaps at the solutions of the scaling 
equation for interface scaled gap. These gaps converge to the interface scaling dimension. 
The results for the interface scaling dimension are shown in Table IIIIl These results are to 
be compared with the known interface exponent in the 0{n) universality class [l^l; which 



are given by the entry p = 1, g = 2 in the Kac table: 

X^ = ^~l. (52) 

which is obtained by the substitution of m = l/{g — 1) in the more common form Xi^2 = 
[(m — 1)^ — l]/2m(m + 1) of the Kac formula. 

Our result for X^ at n = 2 is different from X12 = 1/2 and thus illustrates that the 
n = 2 cubic loop model falls outside the 0(2) universality class. It is to be compared with 
results for the Ashkin- Teller model which are also summarized by Baxter 0. In 

the notation used there, we have the exact result (3^^ = 3/4 for the end point of the Ashkin- 
Teller line at K ^ 00. This exponent may, in our notation, be put equal to Xm/(2 — Xt) 
where Xt = 3/2. This does indeed lead to X^ = 3/8. 

We also calculated the temperature scaled gaps Xt{x,L) at the extrapolated critical 
points. We expect the following behavior of these gaps: 

XtiL)=Xt+puLy- + -- - , (53) 

where p is another unknown constant. Similar fits as before lead to results for the tempera- 
ture scalingdimension that are included in Table ITTl The results agree with the theoretical 
prediction for the 0{n) model 

Xt = --2, (54) 
9 

except for the case n = 2, where our numerical result goes to 3/2 in accordance with the 
exact result for the temperature scaling dimension of the Ashkin- Teller model for 

K ^ 00. 



VI. DISCUSSION 

In general, the discreteness of the n-component cubic model defined by Eq. ^ will enforce 
the existence of a long-range ordered phase at sufficiently low temperatures, also for n > 2. 
However, our transfer-matrix calculations for f = in Eq. ((Tj) did not yield any evidence 
for a phase transition to the ordered phase for n > 2. The absence of a phase transition is 
understandable in terms of the parameters K and M in Eq. (P), because Eq. ((Zj) restricts 
these parameters to -|- M < In 2. This quantity, which specifies the energy difference 
between parallel {si ■ Sj = 1) spins and perpendicular {si ■ Sj = 0) spins, is not sufficient to 
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reach a long-range-ordered phase for n > 2. The maximum value K + M = \n2 along the 
critical line is reached for n = 2 where K = oo. 

For n < 2 we find clear evidence for a phase transition to the ordered phase. In the 
interval 1 < n < 2 we have determined the critical points of the n-component cubic model 
as given by Eqs. (Q) and ((Tj), or by Eq. (jSJ. For the exactly solved cases n = 1 and 
n = 2 we find good agreement with the exact values as given in Table |l] Our results for 
the scaling dimensions, i.e., Xf associated with the temperature, Xh associated with the 
magnetic field, and X^ associated with the introduction of an antiferromagnetic 'seam' in 
the model, agree accurately with the 0(n) universality classes, with the exception of the 
case n = 2 where the model reduces to a special case of the Ashkin- Teller model. For the 
latter case n = 2, our results for the scaling dimensions agree with the exact results for 
the Ashkin- Teller model. The fact that these scaling dimensions are different from those 
for the 0(2) universality class is related with the cubic anisotropy which may be seen as a 
)erturbation of the 0(n) symmetry. This perturbation is irrelevant for n <2 but marginal 
^ for n = 2. This proof of irrelevance for n < 2 applies to small cubic perturbations of the 
isotropy. Our numerical results show that the cubic perturbation remains irrelevant even in 
the extreme anisotropic case described by Eqs. (0) and ((Tj). However, when n approaches 2, 
the exponent i/c associated with the cubic anisotropy field approaches marginality and our 
numerical results thus become less accurate. The cubic anisotropy is truly marginal at n = 2 
and parametrizes the Ashkin- Teller model. For this reason, the cubic crossover phenomena 
are absent forn = 2 and the numerical results are again relatively accurate. The temperature 
scaling dimension of the cubic model and its consistency with 0{n) universality were already 
determined Q] for a range n < 1. The present results extend the range of n and provide 
additional evidence concerning the dimensions Xh and Xm- The numerical accuracy of our 
analysis is such that the 0(n) universality of the cubic model seems reasonably convincing. 
We note that numerical results for the 0{n) model on the square lattice are even more 
accurate; however, for that case the exact critical point is known, and the cubic anisotropy 
field, and thereby the leading corrections to scaling, vanish. 
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